Anisotropic s-wave superconductivity in single crystals CaAlSi from penetration 

depth measurements 
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In- and out-of-plane London penetration depths were measured in single crystals CaAlSi ( Tc = 6.2 
K and 7.3 K) using a tunnel-diode resonator. A full 3-D BCS analysis of the superfluid density is 
consistent with a prolate spheroidal gap, with a weak-coupling BCS value in the a6-plane and 
stronger coupling along the c— axis. The gap anisotropy was found to significantly decrease for 
higher T^ samples. 



Superconductors with AIB2 structure have received in- 
creased attention after the discovery of superconductiv- 
ity at 39 K and especially after identification of two dis- 
tinct gaps in MgB2 0, S 13 ■ K is believed that two gaps 
survive, because of reduced interband scattering due to 
different dimensionality of 2D a and 3D tt bands. Investi- 
gating materials with similar crystal and band structure 
is therefore important for understanding the mechanism 
of superconductivity in this class of hexagonal - layer 
compounds. In this paper we study CaAlSi which has 
been synthesized long ago , but in which superconduc- 
tivity was discovered only recently ^J. Band structure 
calculations show highly hybridized three-dimensional in- 
terlayer and tt* bands [aiLDl- Although, most studies of 
CaAlSi indicate s-wave pairing, deviations from a single 
isotropic gap behavior have been reported |5ll8H9l lldllll |. 
Magnetic measurements indicate a fully developed s-wave 
BCS gap 1^ H^ . Angle-resolved photoemission spec- 
troscopy jlOj revealed the same gap magnitude on the 
two bands with moderate strong coupling value for the 
reduced gap, 2Alk uTr — 4.2. Together with specific 
heat measurements [llj it provided reliable evidence for 
a three-dimensional moderately strong-coupled s-wave 
BCS superconductivity. On the other hand, /iSR stud- 
ies have been interpreted as evidence of either one highly 
anisotropic or two distinct energy gaps Furthermore, 
5-fold and 6-fold stacking sequence of (Al,Si) layers cor- 
responding to two different values of Tc of 6 and ~ 8 
K were found |l3l| . Therefore, an experimental study of 
in- and out-of-plane superfluid density is needed to un- 
derstand anisotropic superconducting gap structure to 
help understanding the mechanism of superconductivity 
in AIB2 type compounds. 

Single crystals of CaAlSi were grown from Ca:Al:Si 
(1:1:1) ignots using a floating zone method as described 
elsewhere Samples have Tc either 6.2 or 7.3 K, which 
is directly related to different stacking sequences 13] . We 



measured two of each types of slab-shaped crystals with 
typical dimensions 0.3 x 0.3 x 0.4 mm"^. The penetration 
depth was measured with an LC tunnel-diode oscillator 
which is sensitive to changes in susceptibility of several 
pico-emu or, equivalently, to changes in London pene- 
tration depth of about 0.3 A for our crystals 17]. The 
quantitative analysis of the frequency shift depends on 
the sample shape and relative orientation of the excita- 
tion field, Hac^ with respect to the principal axes. As- 
suming superconducting crystal with isotropic in-plane 
response determined by the in-plane penetration depth, 
^ab [T] and possible different value of the c— axis pene- 
tration depth, Ac (T), at least two experimental arrange- 
ments are required to extract Xab (T) and Ac (T) sepa- 
rately. In the Hac\\c— axis orientation, superconducting 
currents are generated in the a6— plane, thus the suscep- 
tibility is determined only by Xab (T) and the frequency 
shift. A/ (T) = / (T) ~ /o, is given by 
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where Vq is the effective coil volume, TV is the demagne- 
tization factor, A (T) is the London penetration depth, 
R is the effective planar sample dimension With 
magnetic susceptibility x this equation is just A/ (T) = 
— dTTx (T) A/o where the only sample shape - dependent 
parameter, A/o, is measured directly by pulling the sam- 
ple out of the coil at low temperature. 

Figure ^ shows the in-plane penetration depth deter- 
mined from the frequency shift using Eq.J^I- The zero- 
temperature value was estimated from the fit to the BCS 
formula as described below. This value is not impor- 
tant for the analysis of AAa^ (T), but is needed to esti- 
mate the superfluid density. Muon spin rotation gives 
Xab (0) — 2390 A j^, measurements of the critical fields 
Xab (0) = 2060 A, Ac (0) = 870 A [i2 and Xab (0) = 3140 
A 9] as well as measurements of the reversible magneti- 
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FIG. 1: Low-temperature variation of the in-plane penetra- 
tion depth in two crystals plotted versus reduced temperature 
T/Tc- Red line is a fit to low-temperature s-wave BCS model, 
Eq.0. Inset shows full temperature scale along with the nor- 
mal state resistivities extracted from the skin depth. 
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FIG. 2: c-axis penetration depth obtained from the numerical 
inversion of Eq.||3J for two different samples. Main frame - for 
sample with Tc — 6.2 K and inset shows sample with Tc — 7.3 
K. Solid lines are fits to the low-temperature BCS expression, 
Eq.Q. 



zation Xab (0) = 3100 A 9]. Our conclusions arc not sig- 
nificantly affected by the variation of these values. The 
main purpose of Fig^ is to compare crystals with dif- 
ferent Tc- The inset shows data on a full temperature 
scale. The signal saturates at the level corresponding to 
the normal-state skin depth, thus providing additional 
information - contact-less measurements of the resistiv- 
ity above the transition. We find 45 ^Q-cm and 33 /ifJ-cm 
for 6.2 K and 7.3 K, respectively, which is in an agree- 
ment with direct measurements on single crystals, which 
found 36 /if2-cm on higher Tc sample m ll^. On the 
contrary, when Xab (T) is plotted versus reduced temper- 
ature T/Tc, the curves for two samples coincide (no nor- 
malization was done for the j/— axis). Open symbols in 
Fig. Qshow results for Tc = 7.3 K crystal, whereas closed 
symbols show Tc = 6.2 K material. The data are well fit 
by the standard weak-coupling s-wave BCS model: 



side w the following solution has been obtained |2C 
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where, from AXab {T), we obtained Aab (0) = 1.76A;_bTc, 
- a weak-coupling s-wave BCS superconducting gap. 

In the Hac\\ab orientation shielding currents flow along 
both the a5— plane and the c— axis. The full magnetic 
susceptibility is obtained from the anisotropic London 
equation, which must be solved numerically to extract 
the inter-plane penetration depth Ac (T). For a slab 2b x 
2d X 2w with magnetic field oriented along the longest 
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where kn = tt (rt + 1/2) and 6„ = {{knXab/d)^ + 1^ ■ 

Knowing A^h (T) from independent measurements in the 
i/aci |c orientation and measuring the total frequency shift 
upon extraction of the sample from the coil, A/^""'', 
Eq.© is solved numerically to obtain Ac (T). 

Figure 121 shows the out of plane penetration depth ob- 
tained from Eq.lj^J. The main frame shows the Tc = 6.2 
K sample, whereas the inset shows data for Tc = 7.3 K 
sample. Solid lines are the fits to the low-temperature 
isotropic BCS expression, Eq.j^Jl. The gap amplitude 
obtained from the fits indicates stronger coupling along 
the c— axis (compared to the Aab (T) fits, FigOJ. How- 
ever, as shown below, such fitting significantly underes- 
timates the anisotropy of the superconducting gap. The 
major problem is that Eq.(|2Jl is only valid for an isotropic 
gap that is constant with temperature, restricting its 
range of validity to T /Tc < 0.35. This restriction ex- 
cludes a significant portion of the data. A more com- 
prehensive analysis requires a determination of the nor- 
malized superfluid density, which is obtained from the 
measured change in the penetration depth, AA(T), via 
p(T) = (l + AA(T)/A(0))-^ 

The superfluid density generally depends on the shape 
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of the Fermi surface and the gap anisotropy [21|. For 
CaAlSi we can assume a fairly isotropic Fermi surface |^ 
0, J a superconducting gap isotropic in the ab— plane 
and anisotropic for the out of plane response. Within the 
semiclassical approximation 1215, 
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where z = cos (9) and 9 is the polar angle with 9 = 
along the c— axis. Since there is no general argument for 
the shape of the gap, we choose the spheroidal form: 



T =6.2 K 




A(T,^) = 



v/l-ecos2 (9) 



(6) 



and parameter — oo < e < 1 is related to eccentricity e 
as e = = 1 — c^^ where c is the normalized semi- axis 
along the c— axis, but e can assume both negative and 
positive values. The spheroid is either prolate (e > 0), 
oblate (e < 0) or a sphere (e = 0). The temperature de- 
pendence of the superconducting gap was obtained from 
the anisotropic gap equation. We found that is is well ap- 
proximated by A (T) = A (0) tanh (l.785y^Tc/T - . 

The following data analysis was performed. By mea- 
suring the same sample in two orthogonal orientations 
(along the c— axis and along the a6— plane), both Xab (T) 
and A/c were obtained. The latter contains contributions 
from both Xab {T) and Ac (T) . Equation ^ was then 
used to numerically evaluate Ac (T) . Low-temperature 
BCS fits as well as measurements of the reversible mag- 
netization were used to estimate Xab (0) = 0.31 iim and 
Ac (0) = 0.65 ^m. Fits over the full temperature range 
confirmed the assumed values. Then Eqs. and (O 
were used to fit the data. As a final step both curves, 
Pab{T) and pc(T), were generated from a single set of 
fitting parameters. 

Figure 13 shows data and fitting results for the lower 
Tc = 6.2 K samples. Symbols are the measured data 
points and solid lines are calculated for the ellipsoidal 
gap shown in the inset. The fitting procedure yielded 
the weak-coupling BCS value, 2Aab{0) /ksTc = 3.53, 
in the ab— plane and e = 0.656 corresponding to 
2Ac (0) /ksTc = 6.02 gap maximum along the c— axis. 

Figure ^ shows similar results for the samples with 
Tc = 7.3 K. There is an obvious reduction of the gap 
anisotropy. The best fit to the ellipsoidal gap yields 
e = 0.206 resulting in 2Ac (0) /fc^Tc = 3.98. It should be 
noted that in most previous works only averaged values 



FIG. 3; Two components of the superfluid density for samples 
with Tc = 6.2 K. Symbols show measured data and solid 
lines corresponding full 3D Siting. The dashed line shows 
isotropic s-wave BCS result. The inset shows ellipsoidal gap 
that describes both curves. 
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FIG. 4: Superfluid densities for the higher Tc — 7.3 K samples 
showing apparent reduction of the gap anisotropy copared to 
Fig. El The inset compares cross-section of the gap amplitude 
for the two cases. 



of the superconducting gap could be obtained. From heat 
capacity measurements, 2A(0)//cbTc = 4.07 was ob- 
tained jlli], whereas ARPES yielded 2 A (0) /fc^Tc = 4.2 
[lfij |. The average effective gap can be obtained from 
our results by equating volumes of the spheroidal gap 
and a sphere, Ag// = Aaf, (0) (1 — This gives 
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2Aeff {0) /ksTc = 4.22 for samples with = 6.2 K 
and 2Ae// (0) /ksTc = 3.66 for samples with Tc = 7.3 
K which is in the correct range of reported values and 
our earlier fits using Eq.©. All these values should 
be compared to the weak-coupling isotropic result of 
2A (0) /UbTc = 3.53. 

For all samples studied, we find that the temperature 
dependencies of both in-plane and out-of-plane superfluid 
density are fully consistent with single-gap anisotropic 
s-wave superconductivity. The gap magnitude in the 
a6— plane is close to the weak-coupling BCS value while 
the c— axis values are somewhat larger. Our results sug- 
gest that scattering is not responsible for the difference 
in Tc- Scattering would lead to a suppression of the 
gap anisotropy |23|. The gap with average value A and 
variation (5A on the Fermi surface can only survive if 
hT~^ <C vAi5A, where r is the impurity scattering rate 
[2^ . The values of resistivity are very close for both 
high- and low- Tc samples (45 fjiVl-cm and 33 /zJ7-cm, re- 
spectively), see Fig. the qualitative trend in anisotropy 
is just oppisite. Also, 15% suppression of Tc by non- 
magentic impurities requires very large concentrations. 
This would, indeed, significantly smear the transition, 
which we did not observe. Therefore, all facts point out 
that in CaAlSi gap anisotropy abruptly decreases as Tc 
abruptly increases from ~ 6 to ~ 8 K. A plausible mech- 
anism comes from the analysis of the stacking sequence 
of (Al/Si) hexagonal layers 10. There are two struc- 
tures - 5-fold and 6-fold stacking corresponding to low 
and higher - Tc samples, respectively. Buckling of (Al,Si) 
layers is greatly reduced in a 6-fold structure, which leads 
to the enhancement of the density of states, hence higher 
Tc- Our results suggest that reduced buckling also leads 
to almost isotropic gap function. This may be due sig- 
nificant changes in the phonon spectrum and anisotropy 
of the electron-phonon coupling. 

Measurements of the field dependence of the penetra- 
tion depth in the vortex state also show a difference be- 
tween the two sets of samples, as do measurements taken 
with different field orientations relative to the c— axis. 
Of particular interest is the variation of the second crit- 
ical field and the so-called peak effect HQEIEIEI- 
Tunnel-diode studies of these properties will be reported 
elsewhere. 
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